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Abstract.  The  LS-DYNA  code  was  applied  to  split  Hopkinson  pressure  bar  tests  on  a  material  at  least 
two  orders  of  magnitude  stiffer  in  dilatation  than  in  shear.  Two  constitutive  models  were  applied, 
linear  elasticity  and  a  compressible  form  of  Mooney-Rivlin  elasticity.  The  latter  was  fitted  to  data 
from  ballistic  gelatin.  The  incident  and  transmission  bars  were  aluminum.  The  nominal  strain  rate  was 
2500/s.  Gaps  appeared  at  the  interfaces  between  the  specimen  and  both  bars.  Unloading  of  the 
specimen  and  bars  accompanied  these  gaps.  The  input-velocity  rise  time  was  varied  to  observe  pulse 
shaping  effects.  Mesh  sensitivity  and  contact-parameter  sensitivity  studies  were  performed. 

Keywords:  split  Hopkinson  pressure  bar,  ballistic  gelatin,  LS-DYNA,  interfacial  gaps. 

PACS:  62.20.D-,  83.60.-a,  83.60.Uv. 


INTRODUCTION 

The  LS-DYNA  code  [1]  was  applied  to 
axisymmetric  finite-element  simulations  of  split 
Hopkinson  pressure  bar  (SHPB)  tests  on  a  material 
much  less  stiff  in  shear  than  in  dilatation.  This 
class  has  relevance  to  various  biological  materials 
and  to  their  simulant,  ballistic  gelatin.  The 
measurement  of  high-rate  properties  for  such 
materials  is  challenging  due  to  their  small  wave 
speeds  and  impedances.  A  recent  report  has 
documented  progress  in  extending  the  SHPB 
technique  to  ballistic  gelatin  [2].  The  original 
motivation  for  our  computations  was  to  study  the 
degree  to  which  equilibrium  was  achieved  in  these 
experiments.  Both  linear  elasticity  and  a 
compressible  form  of  Mooney-Rivlin  elasticity 
were  applied  to  the  specimen.  An  unanticipated 
finding  was  the  occurrence  of  separation  gaps 
between  the  specimen  and  the  aluminum  bars  and 
accompanying  unloading  of  the  specimen  and  bars. 


THE  MODEL 
Constitutive  Modeling 

The  aluminum  incident  and  transmission  bars 
were  modeled  with  linear  elasticity.  Aluminum 
was  assigned  the  density  of  2700  kg/m3.  Young’s 
modulus  of  68  GPa,  and  Poisson’s  ratio  of  0.33. 

The  specimen  was  modeled  as  elastic  with  an 
initial  density  of  1000  kg/m3.  First  a  compressible 
generalization  of  Mooney-Rivlin  hyperelasticity 
(“CMR”),  LS-DYNA  model  27  [1,3],  was  applied. 
The  strain  energy  per  volume,  W,  is  related  to  the 
three  invariants  of  the  right  Cauchy-Green  tensor, 
I,  II,  and  III,  by 

W  =  A(l  -  3)  +  B(ll  -  3)  +  C{IIP2  -  1)  +  D(lll  -  I)2  (1) 

Here,  A,B,C ,  and  D  are  constants  related  by 
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_  A(5v  -  2)  +  6(1 1v  -  5) 
“  2(1  -  2v) 


(3) 


vz(t)  = 


(4) 


and  v' is  Poisson’s  ratio. 

In  calculations  applying  CMR,  A  and  B  were  set 
to  12  and  28  kPa,  respectively,  in  order  to  produce 
a  stiffness  intermediate  between  uniaxial-stress 
measurements  reported  in  [2]  at  strain  rates  of  1/s 
and  2500/s.  Constant  v  was  set  to  0.49999  to 
produce  an  initial  bulk  modulus  four  orders  of 
magnitude  greater  than  the  shear  modulus.  In 
calculations  applying  linear  elasticity,  Young’s 
modulus  was  fixed  at  430  kPa  and  Poisson’s  ratio  v 
was  varied. 

Geometry  and  Meshing 

All  simulations  were  2D  axisymmetric.  The 
specimen’s  initial  thickness  and  radius  were  1.45 
and  6.35  mm,  respectively,  the  dimensions  used  in 
[2].  The  bars  had  a  12.8-mm  radius.  The  incident 
and  transmission  bars  had  lengths  of  768  and 
256  mm,  respectively,  or  30  and  10  times  their 
diameter,  respectively. 

All  finite  element  meshes  were  composed  of 
four-node  quadrilaterals.  In  the  bars,  all  elements 
were  1 00-pm  squares.  In  most  of  our  calculations, 
the  specimen’s  elements  were  initially  rectangular, 
having  radial  and  axial  dimensions  of  25  and 
50  pm,  respectively.  Several  calculations  were 
performed  with  different  specimen  meshes  to 
observe  convergence. 

Boundary  Conditions 

Interactions  between  the  specimen  and  the  bars 
was  governed  by  a  penalty-based  contact 
algorithm,  LS-DYNA’s  “2D  Automatic  Surface-to- 
Surface”  [1,3].  The  algorithm’s  key  parameters  are 
SFACT,  a  scale  factor  applied  to  the  penalty  force 
stiffness,  and  VDC,  a  viscous  damping  coefficient 
in  percent  of  critical  damping.  The  defaults  of  1 
and  10,  respectively,  were  used  in  all  calculations 
from  which  results  are  shown.  Parameters  SFACT 
and  VDC  were  varied  in  other  calculations  to 
examine  convergence. 

The  outer  radii  of  the  bars  and  specimen  were 
stress-free.  At  the  input  end  of  the  incident  bar, 
axial  velocity  was  prescribed  in  the  form 


vo~  ;  o<t<tR 
;  t*tR 

Here,  t  is  time.  Rise  time  tR  was  assigned  the 
values  1  and  25  ps  to  explore  effects  of  pulse 
shaping.  The  value  1.81  m/s  was  applied  to  v0  to 
achieve  a  nominal  strain  rate  of  2500/s  across  the 
specimen.  At  the  far  end  of  the  transmission  bar, 
LS-DYNA’s  “2D  Non-Reflecting”  boundary 
condition  [1,3]  was  applied. 


Results  and  Discussion 

Calculations  with  Compressible  Mooney-Rivlin 

Let  AS-ib  denote  the  axial  separation  between 
the  specimen  and  the  incident  bar  at  the  centerline, 
r  =  0.  Let  As_tb  denote  the  axial  separation 
between  the  specimen  and  the  transmission  bar  at 
r  =  0.  Let  /L  denote  the  specimen’s  nominal  axial 
stretch  ratio,  defined  as  the  distance  between  the 
incident  and  transmission  bars  at  r  =  0.  Figure  1 
shows  results  for  these  quantities  vs.  time  for 
tR  =  1  ps.  At  r  =  0,  the  longitudinal  wave  in  the 
incident  bar  arrived  at  the  specimen  interface  at 
126  ps.  By  165  ps,  the  specimen  has  pulled  away 
from  the  incident  bar,  thereby  producing  a  non¬ 
zero  AS-ib-  By  166  ps,  the  specimen  has  begun  to 
lag  behind  the  transmission  bar,  leading  to  a  non¬ 
zero  As_tb.  Three  distinct  intervals  can  be 
discerned  in  Fig.  1.  The  first  interval  occurs 
between  165  and  about  280  ps  and  corresponds  to 
a  A-  range  of  0.97  and  0.68.  In  this  interval  AS_IB 
remains  non-zero  and  attains  a  38-pm  peak  value, 
and  AS-tb  is  intermittently  non-zero  and  attains  a 
peak  of  3  pm.  Both  gaps  remain  closed  in  the  brief 
second  interval.  The  third  interval  extends  from 
about  185  to  325  ps  and  corresponds  to  a  range  of 
0.67  to  0.57  for  Az.  During  the  third  interval,  As_ib 
and  As_tb  are  of  similar  amplitudes  and  attain  peak 
values  of  43  and  39  pm,  respectively. 

Figure  2  compares  As_ib  for  tR  =  1  and  25  ps. 
Pulse  shaping  has  decreased  the  size  of  the  gaps. 
The  peak  value  in  Fig.  2  has  decreased  from  43  to 
3  pm.  However,  even  these  small-amplitude  gaps 


were  accompanied  by  stress  unloading  in  the 
specimen  and  in  the  bars.  This  can  be  seen  in 
Fig.  3,  which  compares  AS.TB  with  axial  stress  <j:z 
in  the  specimen  element  and  in  the  bar  element 
closest  to  the  S-TB  interface  at  r  =  0  for  tR  =  25  ps. 
The  initial  compressive  pulses  that  arrived  in  both 
bars  at  about  150  ps  were  reduced  to  zero  when  the 
gap,  of  only  1.2  pm  amplitude,  formed  at  184  ps. 
Thereafter,  intermittent  gap  closings  were 
accompanied  by  resumptions  in  compressive  stress. 


Such  unloading  in  the  bar  and  specimen  as 
observed  in  Fig.  3  would  invalidate  the  processing 
procedure  usually  applied  to  SF1PB  data  [4].  Axial 
stress  (7-z  is  usually  assumed  to  be  uniform 
throughout  the  specimen  and  is  related  to  the  strain 
8--_  measured  on  the  outer  surface  of  the 
transmission  bar  by 

azz  ~  -jr~E£zz  (5) 

"s 


t  (ns) 

Figure  1.  Gap  amplitudes  at  the  S-IB  and  S-TB 
interfaces  at  r  =  0  and  nominal  axial  stretch  ratio  vs.  time 
for  Compressible  Mooney-Rivlin  and  tR  =  1  ps. 


t  (ps) 

FIGURE  2.  Gap  amplitude  at  the  S-IB  interface  at  r  =  0 
vs.  time  for  Compressible  Mooney-Rivlin  and  tR=  1  and 
25  ps. 


Flere,  E  is  Young’s  modulus  of  the  bars,  and  AB 
and  As  are  the  cross-sectional  areas  of  the  bars  and 
specimen,  respectively.  Flowever,  the  presence  of 
the  S-TB  gap  implies  that  the  contact  area  is  less 
than  the  entire  cross-sectional  area  as  assumed  in 
the  derivation  of  eq.  5.  Also,  the  surface  unloading 
in  Fig.  3  propagates  into  the  interior  of  the 
specimen  and  delays  the  occurrence  of  a  uniform 
state  of  stress. 

The  AS-ib  and  AS.TB  values  in  Figs.  1,  2  and  3  are 
small  relative  to  each  specimen-element’s  initial 
axialxradial  dimensions  of  25x50  pm.  Two  other 
element  sizes,  50x100  and  12.5x12.5  pm,  were 
therefore  applied  to  study  convergence.  Figure  4 
shows  that,  for  tR=  1  ps,  the  AS_IB  results  have 
substantially  converged  during  the  first  interval, 
between  165  and  280  ps.  Mesh  convergence  has 
not  been  demonstrated  for  many  aspects  of  the 
results,  but  the  degree  of  convergence  in  Fig.  4 
supports  the  existence  of  the  gap  phenomenon. 
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Figure  3.  Axial  stresses  on  both  sides  of  the  S-TB 
interface  and  gap  amplitude  at  this  interface,  all  at  r  =  0, 
vs.  time  for  Compressible  Mooney-Rivlin  and  tR  =  25  ps. 


Convergence  with  respect  to  the  contact 
parameters  was  also  studied.  The  simulation  with 
tR  =  1  ps  was  repeated  with  SFACT  equal  to  0.1 
and  10  and  with  VDC  equal  to  20.  Again, 
substantial  convergence  was  demonstrated  for  As_ib 
during  the  first  interval,  but  other  aspects  of  the 
solution  have  not  fully  converged. 

Calculations  with  Linear  Elasticity 

We  applied  linear  elasticity  to  the  specimen. 
Young’s  modulus  was  fixed  at  430  kPa  and 
Poisson’s  ratio  v'was  varied.  The  results  for  AS_IB 
for  tR  =  25  ps  are  shown  in  Fig.  5.  This 
corresponds  to  a  nearly  fixed  shear  modulus  and 
varied  bulk  modulus.  No  gaps  were  detected  at 
either  interface  for  a  v  of  0.48  or  less.  The  largest 
peak  values  (with  respect  to  time)  for  both  AS_IB 
and  As_tb  were  observed  with  a  v  of  0.4990.  Peak 
values  diminished  when  v'was  increased  further. 

Figure  6  contains  a  mesh  plot  for  v=  0.49900 
at  286  pis.  The  values  for  both  AS_IB  and  AS_TB 
exceed  the  specimen  element  dimensions.  This 
observation  further  supports  the  assessment  that  the 
gap  phenomenon  is  not  entirely  an  artifact  of 
meshing. 


CONCLUSIONS 

We  have  not  found  any  explicit  mention  of  the 
gap  phenomenon  in  the  SHPB  literature,  either 
experimental  or  computational.  Each  of  our 
specimen  materials  had  a  shear  modulus  much 
smaller  than  its  bulk  modulus,  an  impedance  much 
less  than  that  of  the  aluminum  bars,  and  was  purely 
elastic.  It  is  not  yet  clear  which  of  these  properties 
was/were  essential  in  producing  the  gaps. 
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Figure  4.  Gap  amplitude  at  the  S-IB  interface  at  r  =  0 
vs.  time  for  Compressible  Mooney-Rivlin,  tR=  1  ps,  and 
three  different  specimen  meshes. 
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Figure  5.  Gap  amplitude  at  the  S-IB  interface  at  r  =  0 
vs.  time  for  tR  =  25  ps  and  linear  elastic  specimens  with  a 
Young’s  modulus  of  430  kPa  and  various  Poisson  ratios. 
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Figure  6.  Mesh  plot  at  286  ps  for  the  linear  elastic 
specimen  with  a  Poisson’s  ratio  of  0.49900,  a  Young’s 
modulus  of  430  kPa  and  for  tR  =  25  ps.  The  centerline  is 
on  the  left  and  the  incident  bar  is  at  the  bottom. 
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1530  3RD  AVE  S 
BIRMINGHAM  AL  35294-4440 

1  PURDUE  UNIVERSITY 

ARNTCS  &  ASTRONAUTICS 
W  CHEN 

701  W  STADIUM  AVE 
W  LAFAYETTE  IN  47907-2045 

1  UNIVERSITY  OF  MISSISSIPPI 
MECHANICAL  ENGRG 

A  RAJENDRAN 
201  CARRIER  HALL 
UNIVERSITY  MS  38677 

2  WAYNE  STATE  UNIVERSITY 
BIOMEDICAL  ENGRG 
CBIR 

K-H  YANG 

818  W  HANCOCK 

DETROIT  MI  48201 

3  UNIVERSITY  OF  CALIFORNIA 
SAN  DIEGO 

MECHL  AND  ARSPC  ENGRG 
M  MEYER 
S  NEMAT-NASSER 
V  NESTERENKO 
9500  GILMAN  DR 
LA  JOLLA  CA  92093-041 1 

1  UNIVERSITY  OF  UTAH 
MECHANICAL  ENGRG 
R  BRANNON 

50  S  CENTRAL  CAMPUS  DR 
SALT  LAKE  CITY  UT  841 12 

1  BROWN  UNIVERSITY 
ENGRG  DIV 

R  CLIFTON 
182  HOPE  ST 
PROVIDENCE  RI  02912 

2  AIR  FORCE  RSRCH  LAB 
RWMW 

J  L  JORDAN 
L  CHHABILDAS 
101  WEST  EGLINBLVD 
EGLIN  AFB  FL  32542 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


4  SANDIA  NATIONAL  LABS 
MR  BAER  MS  0836 
S  SILLING  MS  1322 
ESHERTEL  MS  1185 
PA  TAYLOR  MS  1160 
PO  BOX  5800 

ALBUQUERQUE  NM  97185 

1  SANDIA  NATIONAL  LABS 
T  J  VOGLER 
PO  BOX  969 

LIVERMORE  CA  94551-0969 

7  LOS  ALAMOS  NATIONAL  LAB 

B  CLEMENTS 
E  N  BROWN 
EMAS 
T  MASON 
PRAE 

D  DATTELBAUM 
M  MURPHY 
PO  BOX  1663 
LOS  ALAMOS  NM  87545 

1  LAWRENCE  LIVERMORE 
NATIONAL 
LAB 
D  GOTO 
PO  BOX  808 

LIVERMORE  CA  94551-0808 

1  APPLIED  RSRCH  ASSOC 
SOUTHWEST  DIV 
D  E  GRADY 

4300  SAN  MATEO  BLVD  NE 
ALBUQUERQUE  NM  871 10 


ABERDEEN  PROVING  GROUND 


42  DIR  USARL 
RDRL  CIH  C 
J  CAZAMIAS 
P CHUNG 
D  GROVE 
RDRL  WM 
J  SMITH 
RDRL  WMM  B 
M  BERMAN 
D  HOPKINS 
A  FRYDMAN 
G  GAZONAS 
M  MINNICINO 
B  POWERS 


C  RANDOW 
J  SOUTH 
RDRL  WMM  D 
ECHIN 

C  FOUNTZOULAS 
RDRL  WMS 
T  JONES 
RDRL  WMT 
P  BAKER 
RDRL  WMT  A 
E  HORWATH 
J  HOUSKAMP 
C  KRAUTHAUSER 
S  SCHOENFELD 
RDRL  WMT 
C  HOPPEL 
RDRL  WMT  B 
R  BITTING 
D  FOX 
R  GUPTA 
S  KUKUCK 
RDRL  WMT  C 
T  BJERKE 
M  FERMEN-COKER 
S  SCHRAML 
K  KIMSEY 
D  SCHEFFLER 
RDRL  WMT  D 
S  BILYK 
D  CASEM 
J  CLAYTON 
N  GNIAZDOWSKI 
Y  HUANG 
M  KIM 
R  KRAFT 
BLOVE 

M  RAFTENBERG 
M  SCHEIDLER 
T  WEERASOORI Y A 
K  ZIEGLER 


Intentionally  left  blank. 


